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Abstract 

Cyclic or fluctuating populations experience regular periods of low population density. Genetic bottlenecks during these 
periods could give rise to temporal or spatial genetic differentiation of populations. High levels of movennent among 
increasing populations, however, could ameliorate any differences and could also synchronize the dynamics of 
geographically separated populations. We use microsatellite markers to investigate the genetic differentiation of four 
island and one mainland population of western tent caterpillars, Malacosoma californlcum pluviale, in two periods of peak or 
pre-peak density separated by 8 years. Populations showed high levels of genetic variation and little genetic differentiation 
either temporally between peaks or spatially among sites. Mitochondrial haplotypes were also shared between one island 
population and one mainland population in the two years studied. An isolation-by-distance analysis showed the Fsj values 
of the two geographically closest populations to have the highest level of differentiation in both years. We conclude that 
high levels of dispersal among populations maintain both synchrony of population dynamics and override potential genetic 
differentiation that might occur during population troughs. As far we are aware, this is the first time that genetic similarity 
between temporally separated population outbreaks in insects has been investigated. A review of genetic data for both 
vertebrate and invertebrate species of cyclic animals shows that a lack of spatial genetic differentiation is typical, and may 
result from high levels of dispersal associated with fluctuating dynamics. 
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Introduction 

The geographic structure and genetic differentiation of popu- 
lations are determined by demographic processes and dispersal 
[1]. Of particular interest in this regard are populations that 
regularly fluctuate in density: cyclic species. These include a 
number of species of voles and lemmings, snowshoe hares, grouse 
and some forest Lepidoptera. For these species, population 
densities can fluctuate over orders of magnitude with peak 
densities occurring every 3 (small mammals) to 1 1 (snowshoe 
hares and forest Lepidoptera) years. An interesting phenomenon 
associated with cychc species is that populations often fluctuate 
synchronously over large geographic areas [2,3]. 

Three mechanisms have been proposed as potential causes of 
the synchronous dynamics among geographically separated 
populations of cyclic species: (1) sufficient movement of individuals 
synchronizes populations, (2) widespread weather patterns influ- 
ence the populations, and/ or (3) widespread movement of natural 
enemies synchronize population declines [4]. Dispersal is most 
likely to occur in increasing and peak populations and would 
contribute to the maintenance of synchrony of population 
dynamics and genetic similarity among them. Mechanisms 2 
and 3 would not be as likely to influence genetic structure of 
populations. 



The genetic structure of cyclic populations could be related to 
their dynamics, in several ways. Elton [5] was likely the first to 
propose that cyclic population crashes could reduce genetic 
variation and cause populations to become genetically differenti- 
ated. He also pointed out that genetic variabihty could increase 
with population growth if selection is relaxed during the phase of 
population increase. A third possibility is that density-dependent 
selection or genetic feedback could occur and cause population 
cycles [6] and predictable modifications of the genetic composition 
of populations. Much of the work on genetics and population 
cycles has focused on mammals: voles, lemmings and snowshoe 
hares [7,8,9,10,1 1]. The relationship between genetic change and 
the cyclic population dynamics of northern mammals has recentiy 
been comprehensively reviewed by Noren and Angerbjorn [12]. 
They conclude that genetic variation tends to be high and that 
dispersal is sufficient to overwhelm genetic drift in cyclic 
mammals. Thus dispersal among high-density populations coun- 
teracts the potential influence of inbreeding or selection [13], while 
also increasing the potential for evolution [1] [14]. Some authors, 
however, suggest that intermediate levels of dispersal allow for the 
greatest adaptive divergence among populations ([15] and 
references therein). 

Patterns of genetic variation of species of forest Lepidoptera 
with cyclic population dynamics have received little attention 
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compared to studies of mammals. Social interactions among 
mammals are likely to influence the occurrence and timing of 
dispersal over the population cycle and this is important for 
maintaining high genetic variation and low differentiation [15]. 
For moths, however, movement is likely to be more passive, and 
thus for these animals, patterns of genetic variation among 
populations and years may differ from those of mammals. 

A good example of cyclic dynamics and synchronous outbreaks 
in an insect is that of western tent caterpillars, Malacosoma 
califomicum pluviale (Lepidoptera: Lasiocampidae). Western tent 
caterpillars occur primarily in western North America from British 
Columbia to California. They have one generation a year and feed 
on a variety of deciduous trees. Gregarious larvae form silken tents 
and live in family groups until the last instar. On the Southern 
Gulf Islands in British Columbia, Canada (Figure 1), population 
outbreaks are relatively synchronous and peak population densities 
are reached within a year of each other approximately every 9 
years as shown in Figure 2 ([16,17] and Myers & Cory, 
unpublished data). The conspicuous tents facilitate the monitoring 
of populations even at low density. Tents that represent individual 
families are very sparse in the troughs of the population 
fluctuations and populations go locally extinct in some places 
[2]. Thus genetic botdenecks are likely to occur on a regular basis. 
The gregarious behavior of tent caterpillars probably adds to the 
potential for inbreeding in this species at low densities as well; 
larvae sometimes pupate in the tent and this possibly results in 
siblings mating. This could contribute to the severity of the 
bottleneck and genetic drift. 

Here we use variation in microsatellite markers to test the 
hypothesis that severe bottlenecks and genetic drift during the 
troughs of population density of western tent caterpillars cause 
genetic differentiation among populations or between two 
outbreaks within populations. The alternative hypothesis is that 
dispersal among populations of western tent caterpillars maintains 
synchronous dynamics and genetic similarity across peaks and this 




Figure 1. Map of the Southern Gulf Islands and the lower 
mainland in British Columbia. Islands labeled in bold font were 
sampled for this study. Arrow indicates north. 
doi:1 0.1 371/journal.pone.0096679.g001 



would be supported if populations were not genetically differen- 
tiated. We test these hypotheses by comparing five peak or pre- 
peak populations of western tent caterpillars in 2003 to those eight 
years later in 2011 to determine if different genetic patterns 
emerged following the intervening low density in 2007, if genetic 
variation is maintained, and if island populations are genetically 
differentiated. 

Methods 

Field Collections 

Late instar tent caterpillar larvae and pupae were collected 
during the spring and summer of 2003 and 201 1 when population 
sizes were increasing or had reached peak density (Saturna Island), 
on four islands in southern British Columbia (BC) Canada 
including: Galiano Island (48° 55' 19", -123° 25' 13"), Mandarte 
Island (48° 38' 01", - 123° 17' 14"), Pender Island (48° 45' 07", - 
123° 13' 11") and Saturna Island (48° 46' 32", -123° 07' 55"). 
The fifth population, which we consider to be a mainland 
population, is on Westham Island that is separated from the 
mainland by two arms of the Fraser River (49° 05' 52", — 123° 10' 
19") (Figure 1). Larvae and pupae were brought back to the 
laboratory where the larvae were immediately frozen and pupae 
reared to adult moths, sexed and then frozen at — 20°C or — 80°C 
for later genetic analysis. With the exception of Pender Island, the 
size of the M. c. pluviale population was assessed annually by 
counting the number of tents in the study area or the whole island 
in the case of Mandarte Island, which occupies 7 ha. These 
locations are described in more detail in [16]. Slight variation 
among the population dynamics occurs (Figure 2). The population 
on Saturna Island reached peak density in both 2003 and 2011, 
while populations on Galiano, Mandarte and Pender were at pre- 
peak density those years. The population at Westham was still 
relatively low in 201 1 when the others were high, but it increased 
by 20 1 2 in synchrony with another mainland population that we 
monitored but did not use for genetic analysis. Cypress (Figure 2). 

For the isolation by distance analysis we determined the shortest 
distance between study sites based on their geographic coordinates 
(latitude and longitude) using Google Map. Since these were island 
sites, much of the distance between them is over water. 

DNA Extraction and Microsatellite Analysis 

DNeasy Blood and Tissue kits (QIAGEN, Ontario, Canada) 
were used to extract DNA from larvae and moths. DNA was 
extracted from between 1 1 and 4 1 individuals from each of the 
sites ui 2003 and 2011 for a total of 248 individuals. Only one 
individual was used per tent caterpillar family. Eight microsatel- 
lites that had been proven to be highly variable when character- 
ized in a previous study of M. c. pluviale were selected for the 
analysis of population structure [18]. PCR reactions and 
conditions were as described in [18], with the exception of a 
reduction in the concentration of MgClj from 1.0 mM to 
0.75 mM to reduce stutter for primers A106 and BlOl. 

PCR products were denatured, loaded on a 6% polyacrylamide 
gel, and electrophoresis was performed for 1.5 hours on a LICOR 
4300 automated sequencer (LICOR Inc., NE, USA). Four size 
standards (50-350 bp or 50-700 bp LICOR Inc.), negative 
controls (no DNA) and two reference samples were run on each 
64 well gel to ensure that assigned scores were consistent among 
gels. Gels were scored using SAGA GT version 3.3A for 
microsatellites (LICOR Inc.). 
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Figure 2. Population trends for western tent caterpillars on three southern Gulf Islands and one mainland site. The Pender Island 
population was not counted but was observed to be high in 2003 and 201 1. Vertical lines indicate the years that were sampled for genetic analysis. 
doi:1 0.1 371 /journal.pone.0096679.g002 



Sequencing 

To compare haplotypes and genetic similarity between sites we 
chose two populations for further analysis: Westham Island and 
Galiano Island. Westham Island lies approximately 30 km to the 
northeast of Galiano Island across the Strait of Georgia and the 
two sites are separated by open water. Prevailing westerly winds 
(data from Environment Canada) could assist moth movement 
between these locations and thus similar genetic patterns could 
indicate migration of moths between them. 

We sequenced an intergenic spacer region that had previously 
been found to be highly variable in a number oi Malacosoma species 
(up to 154 bp, Felix Sperling, personal communication). We 
sequenced this region plus part of the COI gene for 12 individuals 
from Galiano Island in 2003, 9 individuals from Galiano Island in 
2011 and 9 from Westham Island in both 2003 and 2011. PGR 
reactions used 10 ng of genomic DNA, 0.2 mM dNTP, (Ix) PGR 
reaction bulfer (Invitrogen), 0.5 units of native Taq (Invitrogen), 
2 mM MgClj (Fermentas, Canada, ON), 5 pmol of forward 
primer (K698: 5' TAG AAT TTA TGG CCT AAA CTT GAG 
GC 3') and reverse primer (K699: 5' AGG AGG ATA AAG AGT 
TGA (C/T)GG 3') (Eurofms MGW Operon, Huntsville, AL, 
USA), and 16.4 |ll HjO for a total reaction volume of 25 |a.l. PGR 
conditions were as follows: 94°G for 2 min, cycled 35 times at 
94°C for 30 s, 45°G for 1 min, and 72°G for 30 s, with a final 
extension of 72°G for 5 min. Products were purified using sodium 
acetate precipitation and sent to Eurofms MGW Operon (AL, 
USA) for sequencing. Bioedit Sequence Alignment Editor [19] was 
used to align sequences and a BLASTn search was performed to 
identify significandy similar sequences. Sequences have been 



deposited in the GenBank database (GenBank Nos. ICF956114— 
KF956152). We used TCS version 1.21 [20] to construct a 
haplotype network based on a combined analysis of the intergenic 
spacer region and partial COI. 

Genetic Diversity 

The number of observed alleles (Na), number of effective alleles 
(Ae), allele frequency, expected heterozygosity (He), and observed 
heterozygosity (Ho) were estimated using GenoDive version 
2.0b22 [21] for each population in 2003 and 201 1. Mean rarefied 
allelic richness (Rx) was calculated in FSTAT version 2.9.3 for 
each population in the two sampling years [22]. 

Population Structure 

FSTAT was used to estimate global Fst values for all 
populations sampled and separately for the two sampling years. 
Pairwise Fst values were calculated in Genodive and significance 
was assessed using 100,000 permutations. Bonferroni correction 
was applied to correct for multiple comparisons. An isolation-by- 
distance analysis was carried out using standardized Fst values 
and log distance between sites using GenAIEx 6.5 [23]. The 
Bayesian clustering program STRUCTURE 2.3.3 [24] was used 
to examine the spatial and temporal structure of M. c. pluviale 
populations. First, simulations were run including all eight 
microsatellite loci for the 248 individuals surveyed from the five 
populations in both 2003 and 2011. Following this, separate 
analyses were performed for each of the survey years. These 
analyses included 124 individuals in 2003 and 124 individuals in 
2011. We followed the recommendations of [25] when selecting 
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parameter settings. For each cluster (A'values) from 1 to 10 we ran 
20 replicate simulations using a burn-in length of 100,000 and 
500,000 Markov chain Monte Carlo (MCMC) repetitions. These 
parameter settings were appropriate given the convergence of all 
summary statistics (a, F, D, and the likelihood). Simulations were 
run using the admixture model with correlated allele frequencies. 
To detect the true number of clusters we calculated the ad hoc 
statistic, ^K [26]. STRUCTURE HARVESTER [27] and 
CLUMPP [28] were used to combine data from replicate runs 
from the same -R value and DISTRUCT [29] was used to visualize 
the results. 

We used TEMPOFs [30] to estimate the variance effective 
population size (NeV) for each of the island populations using 
sampling plan II, where sampled individuals are permanently 
removed from the population prior to reproduction. This method 
requires at least two samples from a population to estimate the 
genetic drift and effective population size between the two dates 
and reduces any bias that might result from small samples sizes or 
skewed allele distribution (at the cost of larger standard deviations). 
The generation time is one year for M. c. pluviale and therefore 
eight generations elapsed between sampling periods. In addition, 
we have compared these data to two other methods which use the 
temporal approach to estimate genetic drift and Ne: the likelihood 
method developed by Wang [31,32] assuming a migration rate of 
zero and the Nei & Taijma [32] estimator Fc, using MLNE [33]. 
We jointly estimated effective population size and migration using 
MLNE [33] for each island population by defining the source 
population based on the pooled allele frequencies for all other 
island populations [34] . 

Results 

A total of 105 alleles were observed for the eight microsateUites 
tested. The number of alleles per locus ranged from 4 (B 121) to 32 
(A106) with an average of 13.12 alleles per locus for the eight loci. 
Seven of the Fis values for each locus for the five islands and two 
years were significant after Bonferroni correction (Critical 
P= 0.05/80 = 0.0006; Table 1). An overaU significant excess of 
heterozygotes is indicated by the Fig value (Fis = 0.07, P<0.001). 
This is due to significant departures from Hardy-Weinberg 
equilibrium observed at three of the loci (B108, B121, and B6) 
(Table 1). These loci do not appear to be influential in other 
analyses and removal of these loci from Est estimates did not 
significantly change the results. 

Sequence Analysis 

Analysis of the intergenic spacer region and partial COl 
revealed 1 1 distinct haplotypes for Galiano and Westham 
populations (Figure 3). Although several of the less common 
haplotypes originated from either Galiano or Westham, there was 
no clear separation of the two sites and the majority of individuals 
were grouped within one dominant haplotype that combined both 
sites and years. 

Genetic Diversity and Population Structure 

The number of alleles observed on each island varied from 6.25 
to 9.75 over both years and after correction for sample size the 
effective number of alleles ranged from 4.54 to 6.33 (Table 2). 
Expected and observed heterozygosity ranged from 0.75 to 0.81 
and 0.64 to 0.71 for populations sampled in 2003 and 2011, 
respectively. None of the diversity measures showed obvious or 
consistent temporal change between the survey years. 

The global Fg-p value for all populations sampled from the five 
islands in 2003 and 2011 was 0.009, which was significantly 
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Figure 3. Network of 11 mitochondrial hapiotypes of western tent caterpillars from Galiano and Westham Islands based on the 
intergenic spacer region and partial sequences of cytochrome oxidase 1. Individuals from Galiano Island 2003 and 201 1 are represented by 
light gray and white and those from Westham Island 2003 and 201 1 are represented by dark gray and black, respectively. The size of each circle is 
proportional to the number of individuals. Lines connecting each haplotype correspond to a single nucleotide change and black lines represent a 
single unobserved nucleotide change between hapiotypes inferred from parsimony analysis. 
doi:1 0.1 371 /journal.pone.0096679.g003 



different from 0 at P<0.001. This indicated a very low level of 
structure among populations. When examined separately the 
pairwise Fg-p values for populations surveyed in 2003 were lower 
than those in 201 1 (Fgr 2003 = 0.007, P<0.05; Fst 201 1 = 0.01 1, 
P<0.001), although none of the pairwise Fg-p values were 
significantly different after Bonferroni correction was applied 
(Table 3). In 2003 five of the Fst values were less than 0.01 and 
three of these involved the Pender population and three the 
Westham population. In 2011 only two pairwise F,st values were 
less than 0.01: one involving Pender and the other involving 
Westham. Isolation by distance analysis showed slight significance 
in 2003 (P=0.04), with the counter intuitive result of Fst 
decreasing (becoming more similar) with increasing distance. No 
significant pattern occurred in 2011 (P= 0.40) (Figure 4). 

In agreement with Fst results, the clustering method imple- 
mented in structure identified litde population structure when all 



populations from 2003 and 2011 were analyzed together. Ln 
Pr(X I A) was highest for the model with K= 1 indicating there is no 
spatial or temporal structure among M. c. pluviale populations 
(Likelihood values: K= 1, -7085, K = 2, -7123, K = 3, -7378). 
Based on the Evanno method, A^STwas highest for K= 3, however 
inspection of the bar plot and individual assignment probabilities 
[q-hat] indicated a low probability of individual membership into 
each cluster for K= 3 and lack of structure for this K value. 
Separate examination of populations surveyed in 2003 (Likelihood 
values: K= 1, -3367, K = 2, -3570, K=3, -3555.6) and 2011 
(Likelihood values: K=l, -3697, K=2, -3931.5, K=3, - 
4094.6) indicated that /S.K was largest for K= 2. However, the 
individuals had a low probability of membership into each of the 
two clusters (cluster 1: 0.50, cluster 2: 0.50) and supported the 
grouping of aU individuals into a single cluster in both survey 
years. 
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Table 2. Genetic diversity estinnates for Malacosoma californicum pluviale populations from five islands in British Columbia, 
Canada. 



Island 


Year 


N3 


A, 


Rt 


Ho 


He 


Galiano 


2003 


6.25 


4.54 


5.98 


0.69 


0.80 




2011 


7.88 


5.31 


6.10 


0.70 


0.76 


Mandarte 


2003 


9.75 


6.33 


6.55 


0.71 


0.81 




2011 


7.5 


4.95 


5.86 


0.70 


0.76 


Pender 


2003 


9.13 


5.16 


6.03 


0.64 


0.75 




2011 


8.13 


5.04 


6.06 


0.65 


0.76 


Saturna 


2003 


9.25 


5.51 


6.27 


0.64 


0.75 




2011 


8.13 


4.96 


5.93 


0.66 


0.77 


Westham 


2003 


7.00 


4.84 


6.13 


0.68 


0.77 




2011 


8.63 


4.99 


6.00 


0.67 


0.76 


Na number of alleles, 


Ae effective number of alleles, 


Rt rarefied allelic richness. 


Ho observed heterozygosity, Hg expected heterogyzosity. 
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Effective Population Size 

Estimates of efTective population size for eacli island population 
based onjorde & Ryman's [30] Fs ranged from 75 to 187, with the 
Pender Island population having the lowest Ne and Galiano Island 
the highest (Table 4). The corresponding upper limits for Galiano 
and Westham Island were infinity indicating that estimates of 
genetic drift could not be differentiated from that arising from 
random sampling error. This occurs when estimates for effective 
population size are low [30]. Fc estimates are within the range of 
values found with other methods, with the likelihood estimate with 
no migration being higher than the others. However, likelihood 
estimates of Ne were low when assuming migration (Ne = 25-50) 
and indicated a high level of migration among all island 
populations (m = 0.9-1.0). This is in agreement with Fst values 
and the STRUCTURE results that indicated a lack of genetic 
differentiation among island populations. 

Discussion 

Estimates of the genetic structure of four island and a mainland 
western tent caterpillar populations indicate low levels of genetic 
differentiation among these populations over scales of tens of 
kilometers. In addition to spatial similarity, population peaks 8 
years apart did not differ. This conclusion is based on the genetic 
similarity of pre-peak populations as indicated by STRUCTURE, 
low Fsx values, the low likelihood estimates of Ne, and the low 
diversity of haplotypes shared between populations. The isolation 
by distance analysis did indicate a weak relationship in 2003 and 
curiously this was the opposite to the anticipated pattern; closer 
populations were more genetically differentiated. There is no 
obvious explanation for this and, given that the result was not 
repeated in 20 11 , further data need to be gathered to assess 
whether this is a consistent pattern. However, it is interesting that 
the two closest populations have almost exactiy the same Est value 
in both 2003 and 2011. 

It is likely that the lack of genetic differentiation is related to 
movement of moths among populations including from the islands 
to the mainland. This suggestion is supported by several lines of 
evidence. First we know that moth dispersal can reestablish "sink" 
populations of tent caterpillars when source populations are high 
[2,35] . In addition, moth movement of a congener, the forest tent 
caterpillar, Malacosoma disstria was observed by Brown [36] to 
occur over 450 km from an outbreak area. Thus movement of tent 



caterpillar moths over tens of kilometers is feasible. The very high 
rate of population increase in the Westham population between 
2011 and 2012 (Figure 2) is compatible with an influx of 
immigrant moths laying eggs there. In addition, as an example 
of immigration we observed a new egg mass at an upper elevation, 
mainland population in July 2012 when the local population was 
still in the larval stages (Myers and Cory, pers obs.). This can only 
be explained by an immigrant moth laying the eggs. Finally the 
typical westerly winds in the evenings and nights when moths are 
flying in the summer in southwestern BC (Environment Canada 
Weather records) would potentially assist movement of moths from 
the island populations towards the mainland. Winds could also 
facilitate movement of moths among islands. The southern Gulf 
Islands are in the rain shadow of the coastal mountains of 
Vancouver Island and the sunnier and dryer conditions there are 
ideal for tent caterpillars. 

In addition to the high levels of genetic similarity among tent 
caterpillar populations, their population dynamics also show 
similar temporal patterns. Tent caterpUlar populations cycle over 
four orders of magnitude and can become locally extinct. The 
resulting population bottienecks could possibly allow for genetic 
drift to occur and influence the ability of populations to adapt and 
evolve. It could also result in genetic differentiation among 
populations if movements were restricted. However, given the 
rapid population increase over three to four generations following 
the trough of population density, genetic diversity should be 
reestablished [37]. 

Estimates of effective population size were less than 200 with the 
use of two of the methods and could not be accurately estimated 
for two of the populations. We found in an earlier study that 
multiple paternity was likely to have occurred in a high proportion 
of families from high-density populations [18]. Waite and Parker 
[38] showed that multiple paternity increases the variance in male 
reproductive success and could influence the estimate of Ne. 

A lack of genetic differentiation appears to be typical of 
outbreak and cyclic species of forest Lepidoptera. Eastern tent 
caterpillars, M. americanum populations were shown with electro- 
phoretic markers to have high genetic homogeneity among widely 
separated populations [39] . Similarly, a genetic analysis of a well- 
studied cyclic forest lepidopteran, the Larch budmoth, Z^iraphera 
diniana, in Switzerland using allozymes showed a low Fsx of 0.002 
among larch feeding populations [40] . The genetic structure of the 
winter moth, Operophtera bmmata, was studied in Belgium using 
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Figure 4. Isolation by distance analysis for western tent 
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allozymes [41] and there was no isolation by distance or genetic 
difiFerentiation on the scale of 10 to 40 km. Winter moth were also 
studied in the Orkney Isles of Scotiand where populations had 
become established in the last 20 to 30 years due to range 
expansion [42]. Genetic differentiation based on AFLP variation 
among populations on three islands was low which suggests high 
gene flow. The population dynamics of this species has been well 
studied in Fennoscandia where cyclic waves of expansion occur 
every 10 years from the east to the west over 3000 km [43]. 
Ballooning of first instar larvae is thought to explain the high 
migration of this species in which females are wingless [42]. A 
review of other studies by Leggett et al. [42] showed that four 
species of Lepidoptera that are known to have high dispersal had 
low Fsx estimates while four species with low dispersal had higher 
Fgx estimates. An example of a species with strong geographical 
structure based on AFLP and mitochondrial markers, and one that 
does not show cyclic dynamics, is pine processionary moth, 
Thaumetopoea pityocampa [44] . The distribution of this pest of pine 
trees has expanded with warming temperatures and pine 
cultivation. Population differentiation is thought to reflect this 
spread and contrasts the results with winter moth [42]. Another 
noncychc, lepidopteran species the cabbage looper, Trichopluda ni, 
migrates aimuaUy to British Columbia from permanent popula- 
tions in southern North America. A high genetic connectivity 
occurred over the migration route, but some genetic differentiation 
occurred among seasonal, green house populations in British 
Columbia [45]. 

Comparisons of the genetic structure of outbreak and non- 
outbreak populations have also been carried out in another insect 
order, the Orthoptera. Outbreak populations of the migratory 
locust, Locusta migratoria, showed much greater genetic homogene- 
ity than did non-outbreaking populations [46]. Similarly the 
plague locust, Chortoicetes terminifera in Australia is characterized by 
outbreak dynamics and showed no genetic structure across the 
country [47]. 

Other non-insect cyclic species also have been shown to lack 
genetic differentiation even though populations periodically 
decline to low density. Burton et al. [7] using eight microsateUite 
markers (the same number as we investigate here) estimated levels 
of genetic variation (13.4 alleles), expected heterozygotes (0.67) 
and Fsx (0.015) for snowshoe hares in the northern Boreal Forest 
of Alaska and Yukon. These values are similar to those we report 
here for cyclic tent caterpillars. They concluded that hares have 
large evolutionary effective populations. Berthier et al. [9] showed 
that cyclic populations of water voles, Awkola tenestris, maintain 
high genetic diversity through migration of increasing populations 
that wipes out the genetic drift among demes that occurs in low- 
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Table 4. Estimates of effective population size (Ne) and migration rate (m) of Malacosoma californicum pluviale populations based 
on temporal samples from surveys in 2003 and 2011. 



Island 


Ne assuming m = 0 






Ne assuming m>0 




Fs 

(Jorde & Ryman 2007) 


Likelihood 
(Wang 2001 ) 


Fc 

(Ne! & Tajima 1981) 


Ne 


m 


Galiano 


187 (77 -=o) 


357 (155-13053) 


176 


28 (16-95) 


0.9 (0.22->l) 


Mandarte 


1 32 (77-470) 


453 (169->10000) 


180 


39 (24-131) 


1.0 (0.18->1) 


Pender 


75 (41-399) 


376 (155->10000) 


198 


50 (28-121) 


1.0 (0.33->1) 


Saturna 


84 (53-91) 


285 (94->10000) 


124 


25 (15-64) 


1.0 (0.26->l) 


Westham 


112 (49-3=) 


256 (78->10000) 


133 


31 (16-332) 


0.96 (0.09->l) 



doi:l 0.1 371 /journal.pone.0096679.t004 



density populations. Similarly Ehrich and Jorde [8] found low 
genetic differentiation and high gene flow using mtDNA markers 
in a study of five species and 72 populations of lemmings in the 
genera of Lemmus and Dicrostorvyx. They suggest that population 
fluctuations possibly destabilize the spatial organization and favor 
the success of migrants. Four-fold variation in mean density of root 
voles (Microtus oeconomus) over six years revealed no change in 
genetic heterogeneity [11]. Rikalainen et al. [10] used 23 
microsateUite loci to study the genetic structure of bank vole, 
Myodes glareolus, populations over eight consecutive years of 
population fluctuation. They found relatively high levels of genetic 
heterogeneity but a lack of genetic differentiation among 
populations at the scale of 10 to 40 km during that time. No 
signature of population botdenecks was apparent. 

Of particular interest with the western tent caterpillars is the 
potential for strong selection associated with the epizootic of a 
nucleopolyhedrovirus that occurs at peak densities which could 
potentially select for different resistant genotypes during sequential 
outbreaks [17]. The variation in disease resistance among 
populations in 2003 is not reflected in the lack of population 
differentiation seen with the microsateUites. It is, however, unlikely 
that neutral markers would show a strong signal of changes in 
other potentially genetically based traits. A similar situation was 
found for cyclic populations of water voles, Arvicola temstris [9] . 
Neutral microsateUite markers indicated low population subdivi- 
sion and low Fgx at high densities. Genes associated with the 
major histocompatibility complex that is related to immune 
function showed evidence of balancing selection and geographic 
diflerentiation among populations at low density. At high density, 
genetic homogenization was greater than could be explained by 
drift and migration alone and therefore could have been associated 
with a general pattern of selection across all populations. 

Genetic variation associated with dispersal of geographically 
separated populations is best studied in metapopulations of the 
GlanviUe fritiUary butterfly, Melitaea cinxia ([48] and citations 
therein). By using functional genomics tools, Wheat et al. [48] 
were able to demonstrate variation between new and old 
populations for loci related to flight and metabolic traits. A more 



genome wide approach such as this will be required to determine if 
ecologically relevant, genetic variation occurs among cyclic 
populations of tent caterpillars. 

We conclude that movement of moths from outbreak popula- 
tions of western tent caterpillars is sufficient to make populations 
genetically similar and to synchronize their dynamics. This pattern 
is typical of mammals with cyclic population dynamics and 
suggests that although opportunities for drift and selection may 
cause populations to differentiate during low density, movement 
among populations at outbreak density eliminates that differenti- 
ation. Further work is required to determine if genetic differen- 
tiation occurs at low density of tent caterpillars, but the lower- 
density population at Westham in 20 1 1 was similar to the other 
populations in that year which indicates that this is not the case. 
Finally, the similar patterns of genetic homogeneity among 
populations of a variety of cyclic species supports the hypothesis 
that dispersal is a common attribute of cyclic animal populations 
and that this maintains genetic similarity and synchronous 
dynamics over wide geographic ranges. 

Supporting Information 

Dataset SI MicrosateUite data for the five Malacosoma califomicum 
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(XLSX) 
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